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Abstract 

The tetrad partial differential equations formulated by Buchman and Bardeen for vacuum 
gravity are shown to be well posed by calculation of the Cartan characters of an associated 
exterior differential system. Gauge specializations are discussed. A Cartan 4-form is found for 
this field theory, together with its intrinsic version the Lagrangian density. 

1 Introduction 

General relativity was discovered and first understood as elegant Riemannian differential ge- 
ometry. For vacuum, this nonlinear theory of gravitation is also a well posed field theory. 
Expressing it most concisely, Buchman and Bardeen jlj have formulated and tested the tetrad 
equations of vacuum general relativity as a first order set for numerical integration. Their depen- 
dent variables are first the connection variables of Estabrook and Wahlquist |2| CHI > a partition 
of the 24 Ricci Rotation CoefRcients Tijk (= —Fi/jj), where «, j, fc = 1...4 into 3-dimensional 
arrays that respect metric signature and local physics and that satisfy directional derivative 
evolution and constraint equations that are symmetric and hyperbolic (the so-called "dyadics" 
Kab, Nabi Wa, Ca, whcrc a,b = 1...3). (An alternate partition of the Tijk is the set of 24 Newman- 
Penrose spinor components.) 

Buchman and Bardeen go on to introduce 16 contravariant coordinate components of the 
orthonormal tetrad field, which we will write ^ , so that partial derivative equations (pde's) can 
be explicitly written. They in fact place both indices to the right, as A^, using Greek indices for 
the range 0,1,2,3. The former notation was used in [2] and [18] and intervening papers to clearly 
distinguish the left index as a label of the tetrad legs (in the terminology of Synge, who put 
such indices in parentheses, it is a "Lorentz" index), while the right is a customary coordinate 
(or tensor) index that is contravariant or covariant, raised or lowered with a metric tensor that 
is a function of position. A Lorentz index is also conveniently raised or lowered, using a (usually 
diagonal) matrix of constants coding the signature of the space. (When there are only Lorentz 
indices and no chance of ambiguity, as in the Tij^ themselves, the indices revert to the right 
side.) In the following we initially take the iX^ (actually, its transposed inverse) as simply a 
matrix of 16 dependent field variables, with i and j labeling rows and columns. Both the tetrad 
orthonormality properties and coordinate covariance properties, the metric interpretations of the 
indices, emerge as a result of the partial differential equations. Buchman and Bardeen's partition 
of these additional fields is a set of dyadic arrays denoted by them , ^b, a, with spatial 



indices that range over 1,2,3. The last two of these are shift and lapse, respectively. Their final 
result is a set of 30 symmetrizeable hyperbolic evolution equations for a total of 40 dependent 
dyadic fields as functions of 4 independent variables x"", t, together with 20 further "constraint" 
equations not involving d/dt. For short we refer to this first order partial differential system as 
the WEBB equations. 

Tetrad and coordinate gauge choices of 10 more evolution equations, together with consistent 
additional constraint equations, remain to be made. Special orthonormal framing due to Nester 

adds six evolution equations for the Tijk and six constraints. The "Lorentz" gauge of van 
Putten and Eardley ^ adds just 6 evolution equations. Harmonic or generalized harmonic 
coordinate specialization can add 4 final evolution equations, for a total of 40, but other 
choices are of course possible. 

Here we use Cartan's theory of Exterior Differential Systems (EDS) to discuss the under- 
lying mathematical structure of the WEBB equations. After setting up an equivalent ideal of 
exterior differential forms in 44 dimensions, we report a calculation of its Cartan characters, a 
table of integers 6, 9 , confirming its Cauchy properties as a well posed field theory. Knowing 
the characters we can count the number of evolution equations, and the number of first class 
constraints, second class constraints, and so on. Cartan characters are also given for different 
gauge choices. We then derive the underlying Cartan 4-form A and its exterior derivative the 
multisymplectic (or Poincare-Cartan) 5-form, and show how the WEBB system results from 
arbitrary variation. They are in a deep sense Euler-Lagrange pde's, but instead of including a 
set of contact 1-forms containing all the partials of the potential fields (as do, e.g., the analogous 
EDS for source- free Maxwell [7] or Yang-Mills [5]) they are built on the torsion 2-forms of Car- 
tan's movable frame bundle, involving at each point a subset of those partials. In the intrinsic 
version of A as a Lagrangian density, the 16 iX^ are potential fields, and the Tijk are intensity 
fields, combinations of partial derivatives of the potentials, and the dynamics, as expressed by 
field equations for Ricci-fiatness, is found by variation. 



2 The Exterior Differential System 

We label 44 coordinates in anticipation of the use to which we will put them as an explicit 
description of 4-dimensional movable frame geometry, i. e. a;* A^- , Fyfe . The basis 1-forms in 
this space are thus dx^ , cPXj, dVijk- In terms of them, we then define sets of 1-forms 6^ and tOij 
(which will shortly be seen as, respectively, tetrad and connection fields) according to 

0^ = ^Xkdx'' 

Lo,j = Vkji ^Xidx^ (1) 

For generality in setting up the EDS we will use a matrix of constants rjij to allow for various met- 
ric signatures. For space-time, and to ensure a hyperbolic system, rjij = rf^ — diag{l, 1, 1, —1). 
In the following all indices except for that on dx^ and the right-hand index on ^Xj will be raised 
or lowered with the t] matrix so that the T]ij need not be explicitly written. Note also that we 
here follow the conventions of and its predecessors in the tetrad components of ojij while 
Buchman and Bardeen [T| label them in the opposite order: TkjiiEW) = Tijk{BB). 

For 0* and ujij to describe a movable frame field and a metric geometry on a submanifold of 
our coordinate space, uJij must be antisymmetric (it is), and the torsion 2-forms must there 
vanish: 

dO' + uj'j A = (2) 

The induced curvature or Riemann 2-forms of the submanifold will be 

c?w,y -f uiik A w J- = Ri-j (3) 

The four Ricci (better, Einstein) 3-forms (which we will want to have vanish) are then A6'fce'-''°' 
and the Ricci "scalar" is the density 4-form Rij A 9k A Oie^^^^ . From these definitions it is easy 
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to calculate 



= (i?^ - uj\ A w*;-) A - Lj'j A {T^ - Lo\ A e'') 

= i?^ A mod r (4) 

and 

dRij ~ Rik A w'^- — LO'^ A Rkj 
d{R^j AOke'^''^) = mod i?y A (5) 

We now can set an ideal of forms (on the 44 dimensional space) whose solutions (submanifolds 
on which the ideal is pulled back to vanish 6,^) will be vacuum metric 4- geometries framed by 
four othonormal 1-forms 6^. This EDS is generated by the torsions, their closure 3- forms and 
the Ricci 3-forms, and we will analyze and interpret it with Cartan theory in the next section: 

{T„i?y A0^i?„ A0fce*-''='} (6) 

These generating forms, it should be emphasized, are all here being given explicitly in terms 
of the 44 coordinates and coordinate bases dx', d^Xj and dTijk ■ This ideal leads to the 50 WEBB 
partial differential equations when one adopts the dx^ as independent variables in a solution, 
inserting dTijk dVijk/dx^ dx^ and d^Xj — > d''Xj/dx^ dx^, after which equating to zero the 
coefficients of the independent 2-forms dx^ Adx^ and 3-forms dx^ Adx^ Adx^ . In a solution the * Aj 
become covariant components of the tetrad frame, and the Tijk its Ricci Rotation Coefficients. 
Buchman and Bardeen more neatly write this system of partial differential equations using the 
equivalent transposed inverse matrix of variables iX^ , contravariant tetrad components related 
to our covariant components by 

= .A'^ d^Xk + 'Xk d^A^ (7) 
Finally, the metric induced on the solutions is 

= m ''A, 'Aj (8) 



3 Cartan Characters, Gauge Freedom and Constraints 

Cartan characters are integers Si, z = 0, 1, 2, ... calculated from the ranks of local linear equations 
for the components of successive vectors Vi,V2,V3, etc., that together span nested solution 
manifolds of an EDS. The Vi are independent up to a limit set by the rapidly growing number 
of equations to be satisfied[Hj|H]. The characters describe the freedom of construction and the 
maximal dimensionality of solutions obtained by successive integrations along a flag of V's 
(thus giving the number of independent variables that can be chosen to write equivalent pde's); 
conversely the characters describe the constraints that arise on successive nested slicings of lower 
dimension, starting from a generic solution ,6j. Simple tests using the characters show whether 
the system is well posed (or causal) , distinguish allowable choices of dependent and independent 
(involutory) variables for writing pde's, and, among the latter isolate the the essential ones from 
those reflecting gauge degrees of freedom. We have calculated the Cartan characters both for the 
WEBB EDS Equation (6) and for versions with various gauge specializations to be given next, 
using Monte Carlo methods developed and programmed in Mathematica by H. D. WahlquistlH]. 

Our notation for the Cartan characters we report places them in a conventional array or 
table N{sq, si, S2, ■■■Sq-i, 0,0,0, ...}g where N is the total number of variables and the Si are 
the successively computed characters. We consider only "well-posed" EDS where for sufficiently 
high i (say q-1) the construction has no more freedom but the sum of the Si has not exceeded 
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N-q. So we set g = — X]?=o understood that there are g-q zeros in the array 

after Sq_i. g is called the genus by Cartan: it is the maximal dimension of an integral manifold 
constructed sequentially as we have described, in practice a generic solution of the associated 
pde system. The zeros indicated that in fact the number of independent variables which can 
be taken is only q — the solution manifold is itself fibered over a q dimensional base, the g-q- 
dimensional fibers express gauge degrees of freedom. Alternatively, we can omit the zeros and 
report, at the last place of the table, q + fiber dimension. Additional gauge choices — more forms 
in the EDS — will result in a lower fiber dimension. In the limit when there is no remaining gauge 
freedom the solution dimension is q=g. For all systems considered here q=4. 

The principal result of this paper is the character array we have calculated for the WEBB 
EDS Equation (6). We find: 44{0,4,12,14}4+10. The EDS passes Cartan's test for being weU 
posed jnj. The number of evolution equations implied is X)i=o ~ '^^^ which agrees with the 10 
gauge degrees of freedom. The total number of pde's will be X^^^q (9 ~ — 50. The number 
of (first class, second class, etc.) constraints is X]i=o ~ * ~ — 20. The number of second 
class constraints is Xli=o ~ * ~ 2)si = 4. 

The 10 fiber dimensions can be understood from the invariance properties of the EDS we 
have set, Eq. (6). An EDS is an ideal of forms, and equivalent sets of generators are related 
by arbitrary invertible algebraic superposition. Six arbitrary tetrad rotations leaving the 77^ 
invariant are allowed due to the homogeneous structure of the Lorentz indices in Eq. (6). Four 
more continuous invariances (or isovectors) are generated by vector fields d/dx^, whose Lie 
derivation acting on the generators vanishes (the do not appear in the EDS except as exterior 
derivatives). Evidently this is the covariance property distinguishing the x'' as coordinates in 
the base. 

4 Gauge Specializations 

The Nester tetrad gauge adds two closed 2-forms to the EDS: 

d(rie^) , d{T,,,eie^^'^') (9) 

The character table computed for this EDS is 44{0,6,14,16}4+4. 36 evolution equations and 26 
constraint equations (6 are second class). 

The Lorentz gauge |l|[Tj adds 6 4-forms to the basic WEBB EDS: 

diT,,tej AOkAeie'^''') (10) 

The character array is 44{0,4,12,20}4+4. 36 evolution equations, 20 constraints (4 second class). 
Since g = 4 and the Lorentz specializations are 4-forms, their closure 5-forms will not change 
the character table, so in fact 6 4-forms like fdx^ Adx^ Adx^ A dec'*, with f an arbitrary function, 
could have been added to them. Buchman and Bardeen pP recognize this in proposing various 
ad hoc right hand sides to the evolution equations for the components of Fyfc that enter Eq. (9) 
and (10), the dyadics and Ga- The left hand sides of these evolution equations are the same 
in both Nester and Lorentz specializations, but the right hand sides in Nester gauge come from 
2-forms, are consistent with the Nester constraint equations, and cannot be arbitrarily changed. 
Harmonic coordinates specialize the ^Xj fields with 4 4-forms: 

dUX'9, ABkABie'^''') (11) 

Note that for conciseness we here have used the contravariant set of tetrad variables iA'*. The 
WEBB character table including these is 44{0,4,12,18}4-|-6. 34 evolution equations, 20 con- 
straints (4 second class). Combined with Nester tetrad constraints one finds 44{0,6,14,20}4 
with no gauge freedom left. 40 evolution equations and 26 constraints (6 second class). Com- 
bined harmonic and Lorentz yields 44{0,4,12,24}4. 40 evolution, 20 constraints (4 second class). 
Since the harmonic conditions are 4-forms the same remark as above applies, that the resulting 
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equations can be generalized with arbitrary functions on the right hand sides of evolution equa- 
tions for the shift and lapse a. Generalized harmonic coordinates have been advocated by 
Pretorius |5|. 

5 Field Theory 

Field theories are usually posited as Lagrange densities, and this can readily be done for the 
WEBB equations. The EDS approach is to formulate a higher rank Cartan form, here a 4- 
form, as a generalization of the Cartan 1-form A = L{q\ (f. t)dt of classical mechanical systems. 
Variation of the Lagrangian L yields canonical equations of motion. Variation of the Cartan 
4-form A is to yield an EDS for the solution geometry and its field equations. In the next 
Section we show how, given A and its EDS, one can finally write a field theoretic, or intrinsic, 
Lagrangian density in terms of partial derivatives. 

Variation of A is achieved by Lie derivation along an arbitrary vector field, and requires that 
the exterior derivative dA (a 5-form) be a sum of terms each of which is at least quadratic in the 
forms generating the EDS. Then the arbitrary variation will indeed vanish on solutions. Such 
an exact 5-form is denoted the multisymplectic or Poincare-Cartan form of the field theory, and 
is a generalization of the symplectic 2-form of classical mechanics. For deep understanding of 
this approach, and bibliography beginning with Caratheodory, Le Page, Dedecker, Rund and 
many others, cf. JO] 

A 5-form can be constructed that is quadratic in the torsion 2-forms and Ricci 3-forms of 
the basic WEBB EDS: 

dK = TiA % A Oue'^'^' (12) 

Each term is the product of a torsion 2-form T/ and a Ricci 3-form Rij A 6ke^^^^ . We have 
denoted it dA anticipating that it is exact, and so a multisymplectic form for the EDS. Indeed 
a quite long calculation, using Equations (4) and (5) and expanding the summations, verifies 
this. c?A can then be integrated by parts to yield the Cartan form: 

A = Rij AOkA 6',e'-^'=' + ^(arbitrary 3-form) (13) 

It is satisfying that we can recognize A as the Ricci scalar, the Hilbert action density for vacuum 
relativity. 

The comment should be made that this factoring of the multisymplectic form dA into a sum 
of products of 2-forms and 3-forms is different from that encountered in most of the literature 
on Cartan forms for field theories [71 |H1 IIUI ITT] The more usual structure is as a sum of prod- 
ucts of 1-forms and 4-forms, with the 1-forms in the EDS being contact forms relating scalar 
potentials and fields and the 4-forms coding the dynamics. Another approach to equations for 
vacuum gravitation that used an EDS for embedding of 4-geometry in flat 10-dimensional met- 
ric geometry ^| led to a Cartan form (a "dimensionally continued" Gauss-Bonnet form |13| ) 
whose exterior derivative is a sum of terms factorable in either of the two ways. The ideal of 
1-forms, closure 2-forms and 4-forms that results from the first factoring {isometric embedding) 
was shown to be well posed; the other factoring, into 2-forms and 3-forms similar to those of 
the present paper, was subsequently considered and is well posed with the same character table 
as that of the basic WEBB EDS jj^. The embedding EDS' are set on a 55 dimensional group 
space (the contraction of the orthogonal group 0(11) to ISO(IO), rotations and translations in 
10 dimensions). Explicit coordinates on the group space are not easy to introduce; also com- 
pared to WEBB there is a larger number of degrees of gauge freedom, viz. 21 (6 for 0(4) tetrad 
rotation and 15 for 0(6) dual tetrad rotation.) The group theoretic embedding EDS may prove 
appropriate for global topological investigations. Alternative well posed field theories are also 
suggested by the embedding formalism [T^ . 
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6 The Intrinsic Lagrangian Density 



To express the variational principle in the more traditional form of a Lagrangian density depen- 
dent on potential fields and their partials, we will begin with a Cartan 4-form A equivalent to 
the above Ricci 4-form for the WEBB EDS, viz. 

A^oji.A uj'^ A 0k A die'^''^ - 2tjy A ujks A 6*" A 9ie'^''' + 2uj,j A A Oie'^''^ (14) 

This form differs from Eq. (13) only by a "surface" term, viz., the exact 4-form d{uoijl\Ok/\Oie^^'^^). 
Its first two terms have been given previously as a "pure grade" Cartan form in the Clifford 
algebra of the second frame bundle over 4-space (where there is no torsion) 15 . Equivalent 
formulations in other notations can be cited, e. g. [Ifil I17| . When we restrict A to a 4-space, 
A Ldx^ A dx^ A dx^ A dx^ , adopting the as independent variables, the resulting Lagrangian 
density is a function of and Fyfc : 

2L = |i?efA,|(rfr^,-r,,fer^"=) (15) 

The first factor is alternatively \DetiX'\~^ . 

To set up Euler-Lagrange equations we need to vary L consistently with the relations between 
and Tijk in the 4-space. Setting T* = 0, i.e., using the EDS to restrict the action in advance 

to a solution 4-space, yields a combination of r^-fe's denoted the "Object of Anholonomity" Q.ijk 

= —flikj in terms of partial derivatives of the *Aj: 

2V,ijk = Tkji — Tjki (16) 
d'Xp/dx" - d'Xs/dx" + 2Ve^k^\J\p = (17) 

The set of 24 Q.ijk are of course equivalent to the Ricci Rotation Coefficients Tijk as they can 
then be inversely solved for: 

Vtsr = ^str ~ ^rts + ^tsr (18) 

Thus in L the 16 ' Aj will play the role of potential fields and the ^ijk or Tijk are intensity fields, 
from Eq.(17) a linear subset ("curls") of partials of the ^Xj. The variations of *Aj and Vijk in 
L must respect Eq.(17), which is 24 of the pde's in the WEBB set (12 evolution equations and 
12 constraint equations for the 'Ap). 16 more WEBB equations (12 evolution equations and 4 
constraints) come from the integrability conditions dTi = 0, the so-called Differential Identities 
satisfied by any tetrad field whether metric or not |18j|19|: 

^r[st.,p\ = '^^\ps^\r\t]q (19) 

We vary L in the customary way using Eq.(17) to express the variations 50, in terms of 
variations of the 'Aj and of the d^Xj/dx^ , permute d/dx^ and 5 and integrate by parts, discard 
surface terms, finding 16 dynamic equations from the vanishing of the coefficients of 5^Xj. They 
are first order in the ^ijk or second order in the 'Aj-: 

'XpdL/d'Xp - 2np^jdL/dnp,m + D„Mldni^^ = o (20) 

where 

Dm =ra X' 3 / dx' (21) 

are directional derivatives in the solution. Inserting the Lagrangian Eq.(15) six of these equations 
turn out to be already included in the set of 16 in Eq. (19). The new ones are the dynamics of 
this field theory, viz. the final 10 WEBB equations (6 evolution, and 4 "energy and momentum" 
constraints), which state in fact the vanishing of the Ricci tensor for vacuum. 
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